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ABSTRACT 

We study the microcanonical density of states and the thermal properties of a 
bosonic string gas starting from a calculation of the Helmholtz free energy in the 
S-representation. By adding more and more strings to the single string system, 
we induce that, for infinite volume, there is no negative specific heat region but 
a transition at a finite value of the energy per string from the low energy regime 
to a region of infinite specific heat at the Hagedorn temperature. Forcing the 
description of this phase in terms of strings gives a picture in which there is 
a very fat string in a sea of low energetic ones. We argue that the necessary 
changing of this description should not change the fact that perturbatively T H is 
a maximum temperature of the system. 
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1 Introduction 


What could be called the Hagedorn problem in String Theory has been a long¬ 
standing issue that lasts for more than twenty five years |IJ. Many of us at some 
moment in the recent past decided to, more or less, forget about it based on the 
belief that non perturbative string theory should have something to say on this 
question. Now that we are able to compute D-instanton effects i i it seems 
very clear that the influence of them on the microcanonical density of states could 
be of some interest. The main question is to know whether there is some sensible 
modification of the picture obtained from the perturbative regime. 

First of all, it would be mandatory to collect all the references about different 
computations of the microcanonical density. There are many of them (cf. for 
example 0P). We focused our interest on the computations we consider as 
more fundamental in the sense that to get the density of states one makes the 
inverse Laplace transform of the canonical partition function. This is actually 
the definition of the microcanonical ensemble when using quantum statistics 0. 
Problems arise soon. The first one is that of choosing a one loop calculation of 
Z(j3) = by taking either what is called the F-representation or the S- 

representation [j7j]. Both are in fact different from one another and only coincide 
in a part of the complex plane once they are analytically continued as functions 
of the temperature. To be more precise, the F-representation does not have a 
unique continuation || to the complex plane and so cannot be used to make 
the inverse Laplace transform. In general, in the F-representation, in which 
modular invariance is manifest, there are many singularities away from the real 
and imaginary axes that do not appear in the S-representation. We have then to 
use the S-representation in which the Helmholtz free energy is computed as the 
sum over the field content of the string. 

Up to our present knowledge, all the calculations in the literature have been 
carried out taking from the very beginning a mostly vague high energy limit 
to obtain an asymptotic form for the density of states which does not contain 
quantum statistic ingredients. Now the question we pose is a little bit naive, 
but to us, will show up as the most important aspect of the problem: how can 
we connect the low energy regime with this high energy behaviour? It seems 
to be clear that a string has two different types of excitations, the massless 
ones which configure what is called the effective theory and all the modes with 
masses multiple of the Planck mass. The massless modes should dominate the low 
energy region and the massive excitations, once the threshold is opened, should 
be responsible for any bizarre behaviour in the high energy part. 

We started thinking about computing D-instanton effects on the density of 
states, but after reviewing the contribution from the perturbative series, we have 
found many things we consider are important enough so as to postpone the 
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original plan. Nevertheless a sketch about the influence of non perturbative 
effects will be depicted in the conclusions. In this paper we shall be devoted to a 
revision of the computation of the string density of states for the closed bosonic 
string using the S-representation. First, we compute the density of states for a 
single string and by repeatedly using the convolution theorem we get the density 
of states for N strings. 


2 Free energy of the bosonic string in 
the S-representation 


The free energy per unit volume for a bosonic string gas in the S-representation, 
also called analog model, is 0 
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where r = T \ + ir 2 . The integral in T \ may be written as: 
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There, we have used the following expansion for the —24 t/l power of the 
Dedekind eta function 


OO 

v ( t )~ 24 = E «ne“. 

n =-1 
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and using an integral expression of the Kd /2 modihed Bessel function, we have 
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Where the fc = —1 term conies from the tachyon and will be neglected be¬ 
cause it produces an infrared divergence; i.e., the tachyon shows its negative 
squared mass when the proper time goes to infinity. The partition function can 
be expanded as 

~ (-W(/3)) n 
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On the other hand, the density of states as a function of the volume and the total 
energy of the system is defined as 
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where C is a Bromwich contour, so the density of states is calculated as the inverse 
Laplace transform of the partition function. If we assume that our system is made 
of many subsystems, that we will interpret afterwards as strings, the density of 
states can also be written as 
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where Qi(Ei) is the density of states of the single string subsystem. We write 
5(x) as a Fourier integral and we obtain 
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Taking the inverse Laplace transform of (|6|), it is easy to see that 
£ _1 {— V/3E(/3)} is the density of states of the single string. This way, the inverse 
of the n th power would be the microcanonical density of the gas with n strings. 
This allows us to relate term by term equations (0) and (0, to give 

co co yn 
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It is worth mentioning that the chemical potential of our system is zero ex¬ 
pressing the fact that the number of subsystems (strings) can change and then 
we have to perform the summation. 

For a certain finite volume, there is always a finite n for which the correspond¬ 
ing term in (|10|) is leading. That would be the most probable number of strings. 
Since in our problem V —> oo, n also goes to infinity, because the last term we 
consider is always the biggest. We will go back to this point later on. 

Usually, when a calculation of the inverse Laplace transform of a product is 
needed, the convolution theorem can be used. Mathematically, it says 

{fi{x)f2(x)} = [ fi(p~t)f 2 (t)dt (11) 

Jo 

in our case: 

1 ME) = 4 (E - t)£h(t)dt (12) 

2! Jo 

that is equivalent, physically, to summing up over every possible distribution of 
the total energy of the system beetwen the two subsystems. Remember that for 
a concrete energy distribution the number of microscopic states is the product of 
the ones of each subsystem. 
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3 Qi(E) at high energies 


A calculation of the microcanonical density of states at high energy can be found, 
for example in ||, obtaining the leading term of the asymptotic series taking only 
into account the ultraviolet singularity of the free energy at the Hagedorn tem¬ 
perature. We have seen that it is possible to obtain the whole series considering 
all the quantum effects, that is, all the Hagedorn type singularities on the real 
axis of the free energy. We write the Helmholtz free energy as in ([!]), and note 
that the part that leads the behaviour of the whole integral at high energy is the 
ultraviolet one, that is: 
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Where the integral of the Dedekind eta function has been approximated in 
the following way. Using the inversion properties of the eta function, it is possible 
to write 



The same Fourier expansion of the —24 th power we used in (||) yields 
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Now, we approximate the series by its first term, defined by n — m — — 1, 
because it is the only divergent one as r 2 —> 0. Expanding the integrand around 
t 1 = 0, which is the region which encloses most of the area in the integral: 


b (t 2 ) ~ ^ 2 rf /2 ( 16 ) 

which is a better approximation as e approaches 0. To finally evaluate the re¬ 
sulting integral in r 2 it is useful to change r 2 by l/r 2 , and we also eliminate 
the tachyon term. Using the integral representation of the incomplete Gamma 
function, the following expression for the single string partition function is found 
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where f3 n = When n = 1 we have the inverse of the Hagedorn temperature 

that we will call fin from now on. Continuing each term of this series to the 
complex p plane it is easy to see that it has branch points at every p n , but every 
term of the series is well defined with isolated singularities at the branch points 
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and the corresponding branch cuts. This allows us to invert Laplace term by 
term looking at each term as an independent inverse Laplace transfom problem, 
what guarantees the convergence of the resulting series, which is 


where 
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and 'H(x) is the Heaviside function. This is a finite series in the n index and an 
asymptotic one in m. The leading term is the well known Hagedorn one, that 
accounts for the Maxwell-Boltzmann statistics. The series does not fundamen¬ 
tally change the behaviour of this term, but it corrects the apparent divergence 
that appeared in previous calculations at low energy because the step functions 
eliminate its contribution there. After all, at high energy Bose and Maxwell- 
Boltzman statistics are expected to be equivalent. Nonetheless this density of 
states of the single string cannot be used to obtain the multiple string through 
convolution because it does not have the information of the low energy range 
which is completely contained in the infrared part neglected in the above cal¬ 
culation by introducing the cut off e. It is important to notice that although 
the behaviour of a function in the complex plane is nearly completely defined by 
its singularities, the regular parts do play an important role when we Laplace 
transform. As an example, in the above calculation the step functions that came 
from the analytical parts guaranteed the convergence of the series in the n index. 


4 Complete single string density of states 

In this section we will obtain the complete single string density of states, and 
compare it with obtained in (|TE|). We will start from equation (|5|). To 

perform the inverse Laplace transform we use the formulas in Jll| to get 
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Remebering that C 1 {/3} equals minus the derivative of the Dirac Delta func¬ 
tion we obtain 
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The most innovative thing is that in our calculation it is explicit the relativistic 
nature of the theory. The Heaviside step functions represent the abrupt appearing 
of the mass levels because while there is not enough energy to create massive 
holds, those states are not accesible. The Hagedorn singularity, explicit in the 
other calculation, is here hidden in the enormous growth of the a\ coefficients, 
that represent the degeneracy of the mass levels. This expression for the density 
of states is valid over the whole range of values of the total energy of the system; 
for example, it possible with it to study the massless quantum held regime of the 
theory and the transition to the ’’Hagedorn phase”. In fact, for energies lower 
than the Planck mass the series has only the k = 0 term so that the system 
behaves exactly as a collection of massless quantum helds. With the density of 
states in eq.(|2T|) the multiple string microcanonical density can be gotten using 
convolution theorem so as to numerically compute it. Notwithstanding, in the 
case of the effective theory, it is easy to obtain a simple expression. The single 
particle density of states is 

Sil ,«//(£) = p 5 P T 7 5 C(26)(24 y-E™ (22) 

where £(;?) is Riemman’s Zeta function. Note that a q = 24 2 represents the degen¬ 
eracy of the massless level. Twenty four is the number of transverse dimensions. 
If we make N convolutions: 

1 / v \ N 

= U25iV)iV! (2W C(26)(24)2r(25) ) E25N ~' (23) 

which is the micraconical density for N massless particles in 25 spatial dimensions. 
We can use this expression as a model to understand the relation between the 
volume and the number of subsystems -in this case, particles. For a given volume 
we can calculate, in terms of N, the most degenerate configuration of the gas, 
which is the most probable. We obtain the most probable value of the number of 
subsystems by evaluating N which is the most relevant term of the series obtained 
summing up all the possible iV’s: 

N oc u 25 V (24) 

where u is the energy per particle. If we take the volume to infinity, the number 
N goes also to infinity, but the density of particles per unit volume remains finite 
as long as the energy per particle is finite. 
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Figure 1: Degeneracy of the two string system, in the broken-equipartition 
regime. The total energy of the system is . 


5 The problem of energy equipartition 


To study the most probable distribution of the total energy of the system among 
the strings we will write the density of states for a given energy distribution. 
Let us imagine a system with two strings that share a total energy E, when one 
of them has an energy £, the microcanonical density is 


u 1 (E,£) = ^Sl l (E-£)(l l (£) 


(25) 


which is the integrand and combinatorial factor in (|12|). The probability would 
be 

^i(E — £)Qi(£) 


V(E,S) = 


V2{E) 


(26) 


Di (E) is the single string density of states that we have already calculated. 

Loocking at this distribution in terms of £ for a given E we see that the behaviour 

is very different whether the total energy is bigger or smaller than E 0 = 2uq ~ 
2.8 
\foJ 

If the total energy is smaller than the mentioned value, the most probable 
distribution is £ — y and then we have equipartition. The strings will most 
probably be in the massless level or in the first Planck mass excitation because 
for these two levels the degeneracy is not high enough to destabilize the system 
for energies lower than E$. But when E > Eq one of the strings tends to absorbe 
more energy than the other, (see Fig.[l|). In this last regime equipartition is 
broken. Another way to observe this behaviour is by looking at the density 
of states of the two string gas obtained through convolution, and comparing it 
with the single string one. As the energy grows both functions give the same 
termodynamics in the sense that for the same big total energy both densities 
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Figure 2: Degeneracy of the three string system with a total energy 


of states produce the same thermal properties, seeQ Fig.||. In the next section 
we will find T as a function of the energy. The same behaviour can be seen 
for a three string gas, as in Fig.^. One could wonder whether the fat (energetic) 
string would eventually capture all the energy while the others disappear. To this 
purpose we can make a simple calculation for the two-string example. We start 
from a concrete configuration of the gas with a total energy E. We assume that 
one of the strings is energetic enough to approximate well its density of states by 
the leading behaviour at high energies, and the other will be supposed to be in 
the massless range. The degeneracy of this configuration would be 

u(E, £) oc £ 2 \E - £)-27/2 e ME-E) (27) 


where £ is the energy of the light string. As a function of £, u(E, £) will increase 
if 


24 27 
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> 0 


(28) 


and decrease otherwise. For a given value of £ if uj(E, £) grows (lessens) the most 
probable value of £ would be bigger (smaller) than the considered one. Note that 
if £ < then the condition will always be satisfied no matter how energetic the 
other string is. This means there will always be a small string of at least that 
critical energy. This simple model allows us to induce that in a string gas there 
will exist, as the total energy becomes bigger, a very energetic string and a low 
energetic string sea. In the next section we will make another explanation of this 
fact for a n-string gas using its thermal properties. 

1 In the plots, a! = 1 has been taken and the temperature has been scaled by 47r so Th = 1. 
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Figure 3: Temperature of the one and two-string gases with the same total energy. 
Both systems have E = 2 u, where u is the energy per string in the two-string 
case. 


6 Thermodynamics 


We start from some basic definitions of thermodynamical magnitudes in the mi- 
crocanonical ensemble: 
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The expressions above for the temperature clarify the problem of the infinite 
volume. Writing Q(E) as a series in n , the temperature becomes 
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Then the only term one has to consider in order to calculate the thermal prop¬ 
erties is the last one, that is the term with the biggest n = N. In other words, 
we are taking V to infinity irrespective of the value of E. 

In general the specific heat is a linear function of the number of subsystems 
and so it is sometimes useful to define a specific heat per subsystem, as follows 


C v l {u) 
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We have computed numerically the temperature for the single and the two- 
string gases. The results are shown in Fig.|3] and Fig.[|. 

Some important features come up. First of all, one sees that both systems 
go from the low energy regime to the high energy one through a maximum. The 
explanation for the single string is clear; when the system is in the low energy 
regime everything happens like in a particle gas, that is, the specific heat is finite, 
positive and constant. However, when the amount of energy is enough for the 
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Figure 4: Temperature of the single and two string gases with the same energy 
per string. 


string to reach the Planck mass level, this energy is in part used to create the 
mass and not to increase the kinetic energy so the temperature tends finally to 
decrease. In a system with only two levels the temperature would increase again 
because after creating the mass the rest of the energy is employed in accelerating 
the particles, and so elevating the temperature. In our system there is an infinite 
number of highly degenerated mass levels, each one more than the previous, such 
that the energy will always be used to create ever more massive fields. That is 
why the temperature gets lower. 

To calculate the high energy limit of the temperature, it is useful to approx¬ 
imate the density of states by the leading term in ([18]). It gives the well-known 
expression for the temperature: 


T 


E 

0hE - f 


(32) 


which tends to T h when E —> oo. 

In the two-string gas the behaviour is practically the same except for the fact 
that the maximum occurs at a lower temperature. This comes from the growth 
of the number of energy levels that are accesible to each string when it borrows 
some energy from the other. The deviation from equipartition is responsible for 
the cooling of the gas. The maximum appears approximately at a total energy 
2«o, that is, when each string has the critical energy assuming equipartition. In 
both figures, we can see the two different regimes of the gas. When the energy 
is large, the two-string gas behaves nearly as a single string with all the energy, 
just as the energy distributions had predicted. 

It is possible to generalize this high energy behaviour for a multiple string 
gas using the degeneracy of a configuration similar to the one in eq.(^). Let us 
suppose we have N strings in the low energy regime and a fat one sharing a total 
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energy E: 


u>(E, u) — Qn(Nu)Qi(E — Nu). (33) 

Then taking the derivative with respect to u, which is now the density of 
energy per string, it is easy to see that u>(E, u ) will increase with u only if Tf at > 
T sea . Since when the energy of the energetic string grows, the value of Tf at tends 
to the Hagedorn temperature Th, then in this situation the temperature of the 
other strings will also be Th- It seems that the string sea plays the role of a great 
heat bath in which the fat string is contained. 

We can put together all the information we have: the single and two-string 
gases density of states, the energy distribution of any number of strings and other 
physical hints that allow us to present a pretty realistic model. 

Let us imagine we have an infinite universe filled with a very rarefied and low 
energetic string gas. Let us suppose we are able to introduce more and more 
energy at the same time that we observe the system evolution. It would begin 
behaving like a particle gas with a positive specific heat per particle. As the 
density of energy per string grows, the strings are able to occupy more and more 
massive states, increasing the specific heat, as we have already explained for the 
two-string gas. At a certain value of the energy, the temperature reaches the 
Hagedorn one. Until that point, equipartition applies, but if we tried to increase 
slightly the energy of the gas, the probable distribution would change in such a 
way that a single string would come out from the gas and absorb all the energy 
added. In other words, the most probable number of strings absorbing energy is 
one. 

It is very important to note that since we are considering an infinite volume 
universe, the number of strings and the total energy of the system are also infinite. 
Hence, it is only possible to deal with densities per unit volume or per subsystem. 
Both densities are simply related at least in the particle case, as we have seen 
before. For energies higher than the critical one, the density of strings is always 
the same as no energy is absorbed or emitted by the string sea. 

It is clear that the fat string appears with an already infinite amount of energy, 
and it can only exist being infinitely energetic. Therefore its temperature will 
always be Hagedorn. We find no negative specific heat for any value of the energy. 
We do, however, find a divergent specific heat for energies per string bigger than 
a certain critical point (approximately Uq). This is not a surprise because when a 
finite maximum temperature of a system is expected, the specific heat diverges. 
For example, the only difference with the open string case is that this divergence 
in Cy appears at a finite value of u. One could wonder whether this divergence 
appears abruptly or in a continuous way. We have seen that when one makes one 
convolution, the maximum temperature gets lower. The corresponding maxima 
are smooth and no discontinuity in T is observed. That is why we think that 
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Figure 5: Tentative plot of temperature versus energy per string 

the natural extrapolation for a large number of convolutions is that the system 
will keep this behaviour, and the maximum, that we already know is Hagedorn, 
will be reached smoothly. We present a tentative plot of the qualitative thermal 
behaviour of the system in Fig.^|. 

7 Conclusions 


In this paper we have revisited the thermal behaviour of a bosonic string gas in the 
microcanonical ensemble. We have introduced some mathematical tools which 
seem to be useful to make the inverse Laplace transorm. The role of the analytical 
parts have manifested as important to grant the convergence of the series which 
allowed us to account for the quantum statistics effetcs in the single string density 


of states at high energies. However, the physical range of application of (|18D is not 
known because it does not contain the infrared part of the free energy associated 
in the (3 complex plane with singularities on the imaginary axis. On the other 
hand this expression is completely useless in order to calculate the multiple string 
density. This is due to the need of the convolution theorem, that we see as the 
most physical way to compute the multiple string case. This enforces us to 
try to compute a valid expression for the density of states for the whole energy 
range. We succeeded making that for the single string and, by application of 
the convolution theorem, we observe a behaviour from which we induce that for 
infinite volume, summing up over an infinite number of subsystems, no negative 
specific heat phase appears. On the contrary, at high energies, we have a system 
which absorbs any amount of energy without increasing its temperature, so the 
limiting temperature is actually a maximum one for the universe. What 
seems to be clear to us is that the description of this phase as one with, not only 
a fat string, but also with a sea of a large number of low-energy strings is not 
adequate to the degrees of freedom we have at hand. It is only the result of trying 
to enforce a description of these degrees of freedom in terms of strings when no 
equilibrium is posible among them as subsystems. To be more concrete, we do 
not believe an alternative description of these degrees of freedom should change 
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the fact that, perturbatively at least, Hagedorn is the maximum temperature of 
the system. 

The study we have presented can be easily extended to Superstrings. The 
next refinement of the thermodynamics would be considering non-perturbative 
effects as D-Instantons. A calculation of the influence of these effects on the 
one loop free energy has been done in J2|, where the authors have also obtained 
an expression for the high energy density of states for this case. This type of 
calculation is along the lines of the works in ff|] in which no infinite volume limit 
is actually taken, but what they do amounts to taking a high energy limit for each 

in (|10| ) then substituting the energy dependence for that of the single string 
in such a way that finally they approximate Q(E) by the asymptotic expression 
for the single string density times the exponential of V. So no limit in V is in 
fact needed. In this way no phase transition can then be captured because no 
critical limit is taken. Obviously this procedure is completely different from ours 
ab initio. 

We believe that to include D-Instantons effects in our picture, it is necessary 
to compute once more the complete density of states for the single string with 
this corrections, and then, by convolution, obtain the multiple string gas. 


Acknowledgements 

We are greatful to J. L. F. Barbon and M.A. Vazquez Mozo for discussions. 
Thanks to E. Alvarez for reading the manuscript and to J.M. Noriega Antuna 
for helping us to find an indispensable book. Jesus and Marco are in debt with 
their parents for the economic support. 


13 


References 

[1] R.Hagedorn Nuovo Cimento Supp.3 147 (1965) 

S.Frautschi Phys. Rev. D3 2821 (1971) 

R.Carlitz Phys. Rev. D5 3231 (1972) 

N.Cabbibo and G.Parisi Phys. Rev. B59 67 (1975) 

[2] J.Barbon and M.A. Vazquez Mozo. CERN-TH/96-361 and IASSNS-96/127 

[3] M.Green Phys. Lett. B354 271 (1995), |icp-th/95041~08| . 

[4] N.Deo, S.Jain and C.-I.Tan Phys. Rev. D45 3641 (1992) 

M.Axenides, S.D.Ellis and C.Kounnas Phys. Rev. D37 2964 (1988) 
M.McGuigan Phys. Rev. D38 552 (1988) 

[5] E.Alvarez Nucl.Phys. B269 596 (1986) 

[6] G.Horowitz Comm. Math. Phys. 89 117 (1983) 

[7] B. McClain B. Roth Comm. Math. Phys 111 538 (1987) 

D. Kutasov and N.Seiberg Nucl. Phys. B358 600 (1991) 

K.H.O’Brien and C.-I.Tan Phys. Rev. D36 1184 (1987) 

E. Alvarez and M.A.R. Osorio Nucl. Phys. B304 327 (1988) 

M.A.R. Osorio Int. J. Mod. Phys. A7 4275 (1992) 

[8] M.A.R Osorio and M.A. Vazquez Mozo Phys. Lett. B280 21 (1992), Phys. 
Rev. D47 3411 (1993) 

[9] E.Alvarez and M.A.R.Osorio Phys. Rev. D36 1175 (1987) 

[10] R.H.Brandenberger and C.Vafa Nucl. Phys. B316 391 (1988) 

[11] A.Erdelyi (Editor) Tables of integral transforms. Vol I (McGraw-Hill Book 
Company, Inc. New York 1954) 


14 




